Models of universe with a polytropic equation of state 
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We construct models of universe with a generalized equation of state p = (cup + kp 
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having a linear component and a polytropic component. The linear equation of state p = ape 2 
describes radiation (a = 1/3), pressureless matter (a = 0), stiff matter (a = 1), and vacuum 
energy (a = —1). The polytropic equation of state p = kp 1+1 ' n c 2 may be due to Bose-Einstein 
condensates with repulsive (k > 0) or attractive (k < 0) self-interaction, or have another origin. In 
this paper, we consider positive indices n > 0. In that case, the polytropic component dominates in 
the early universe where the density is high. For a = 1/3, n = 1 and k — — 4/(3pp), we obtain a 
model of early universe describing the transition from a pre-radiation era to the radiation era. The 
universe exists at any time in the past and there is no singularity. However, for t < 0, its size is less 
than the Planck length lp = 1.6210 _35 m. In this model, the universe undergoes an inflationary 
expansion with the Planck density pp — 5.16 10" g/m 3 that brings it to a size ai = 2.61 10 _6 m at 
ti = 1.25 10~ 42 s (about 20 Planck times tp). For a = 1/3, n = 1 and k = 4/(3pp), we obtain a 
model of early universe with a new form of primordial singularity: The universe starts at t = with 
an infinite density and a finite radius a = oi. Actually, this universe becomes physical at a time 
U = 8.32 1 -45 s from which the velocity of sound is less than the speed of light. When a>fli, the 
universe evolves like in the standard model. We describe the transition from the pre-radiation era 
to the radiation era by analogy with a second order phase transition where the Planck constant fi 
plays the role of finite size effects (the standard Big Bang theory is recovered for Ti = 0). 



I. INTRODUCTION 

The discovery of the expansion of the universe and the 
concept of a singularity from which the universe emerged 
(Big Bang) has an interesting history P, . In this intro- 
duction, we provide a brief review of the early develop- 
ment of modern cosmology, describing very elementary 
(yet fundamental) cosmological models that will play an 
important role in our study. 

In 1917, Einstein 3| applied the equations of general 
relativity to cosmology, assuming that the universe is ho- 
mogeneous and isotropic. On the ground of cultural and 
philosophical beliefs, he considered a static universe. In 
order to accommodate for such a solution, he had to in- 
troduce a cosmological constant A in the equations of 
general relativitjO- He obtained a model of the uni- 
verse with a positive curvature k = +1 that is finite 
(yet unbounded) with a "radius" = c/VA, a den- 
sity p E — A/47rG, and a mass M E = 7rc 3 /2Gv / A. The 
same year, de Sitter 0, [f| discovered another static so- 
lution of the Einstein equations (according to the review 



1 As a preamble to his paper, Einstein Q considered the New- 
tonian approximation and replaced the Poisson equation by an 
equation of the form A<E> — (A/c 2 )'l> = AnGp, so as to obtain a 
static solution <t> = —4nGpc 2 /A. In Einstein's belief, the Newto- 
nian effect of the cosmological constant was to shield the gravita- 
tional interaction on a distance c/y/A. This is actually incorrect. 
Lemaitre [|[ was the first to understand that the cosmological 
constant can be interpreted as a force of "cosmic repulsion". In 
the Newtonian context, the modified Poisson equation including 
the cosmological constant is A<E> = AirGp — A, leading to the 
static solution p = A/AttG. 



of Robertson this solution was suggested to de Sit- 
ter by Ehrenfest). This universe is empty (p = 0) and 
has a radius ads = Cy/3/A. In this model, the stars and 
nebulae are just "test particles" that do not curve the 
universe. De Sitter managed to predict that the light 
emitted by these sources should be redshifted because 
(i) they move away from the observer due to the effect 
of the cosmological constant, and (ii) there is a slow- 
ing down of time with increasing distance from the ob- 
server (de Sitter effect). He mentioned that some obser- 
vations (presumably those made by Slipher since 1912) 
seem to confirm that claim. He also predicted a linear, 
or quadratic, relationship between the redshift and the 
distance. However, his interpretation of the redshifts re- 
mained enigmatic for a long time. Therefore, the early 
cosmological debate revolved around two deficient mod- 
els: The model of Einstein which was unable to explain 
redshifts, and the model of de Sitter which contained no 
matter and was rather mysterious. Both were built on 
the assumption that the universe is static. 

In 1922 and 1924, Friedmann [H [i| discovered non- 
static cosmological solutions of the Einstein equations de- 
scribing the temporal variation of space. He considered 
the case of spherical (positive curvature k = +1) and 
hyperbolic (negative curvature k = — 1) universes. He 
described expanding and contracting solutions, as well 
as cyclic solutions, and foresaw the possible beginning of 
the universe in a singularity. He obtained for the first 
time an estimate of the age of the universe of the order 
of 10 billion years. At that time, Einstein did not believe 
in the significance of these non stationary solutions. He 
first argued that the Friedmann solutions were not com- 
patible with the field equations of general relativity [l(| , 
then retracted his statement II ll but remained convinced 
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that the Fricdmann solutions are not physically relevant 
P, [2[ . Apart from Einstein, who criticized them, the pa- 
pers of Friedmann remained unnoticed by the physical 
community until 1930. 

In 1922, Lanczos [l2| found a non-stationary form of 
de Sitter's solution in which the scale factor increases 
exponentially rapidly. He had the key to an expanding 
universe but did not recognize its physical significance. 
In 1923, Weyl [l3| published a paper in which he tried to 
find his own way to de Sitter's redshift. He noted that 
the particles in de Sitter's universe move away from each 
other with a velocity proportional to their distance, but 
did not provide an explicit system of coordinates to make 
his claim more precise. In 1925, Lemaitre [TiJ found a 
system of coordinates in which de Sitter's model is non- 
static and leads to a linear relationship between the ve- 
locity (redshift) and the distance. This model explains 
the redshifts of distant objects in an unambiguous way. 
This model is Euclidean k = (which was problematic 
for Lemaitre) and the scale factor increases exponentially 
rapidly with time due to the cosmological constant. 

In 1927, Lemaitre [l5| rediscovered independently from 
Friedmann the non-stationary cosmological solutions of 
the Einstein equations. In addition, having spent two 
years in Harvard and at the MIT, he was aware of cos- 
mological observations and related the expansion of space 
predicted by the theory of general relativity to the reces- 
sion motion of galaxies observed by Slipher, Hubble and 
Strdmbcrg. He interpreted the redshift of the nebulae as 
a consequence of the expansion of the universe, predicted 
a linear relationship between the recession velocity and 
the distance (the nowdays called "Hubble relationship" ) 
and, using available data [l6j . obtained an estimate of 
the "Hubble constant" H = a/ a = v/r two years be- 
fore Hubble \YTln. The work of Lemaitre, like the one of 



his predecessor Friedmann, was still not appreciated by 
Einstein who found it "abominable" Q], 3 ■ 

In 1928 and 1929, Robertson [H reinterpreted de 
Sitter's solution dynamically. However, like Lanczos [l2[ 
and Weyl [l~3|, he did not "see" the expanding universe. 



2 As pointed out by Luminet pj and Nussbaumer & Bieri Q, the 
Hubble law and the discovery of the expansion of the universe 
should be attributed in part to Lemaitre, especially if we consider 
that Hubble did not relate his observations to the expansion of 
space in the theory of general relativity (and was later reluctant 
in accepting this idea), while this was Lemaitre's main motiva- 
tion. Indeed, Hubble interpreted the redshifts in terms of de 
Sitter's effect in a static universe, while Lemaitre was the first to 
connect them to the expansion of the universe. Unfortunately, 
the paper of Lemaitre was totally ignored at that time because 
it was published in French in a rather inaccessible Belg ian jour- 
nal. His paper was translated in English in 1931 llSl . but for 
unknown reasons, the important paragraph of the paper where 
Lemaitre gives the relation of proportionality between the reces- 
sion velocity and the distance, and extracts the Hubble constant, 
is condensed into a single sentence. Because of this omission, 
Lemaitre is not recognized on the same footing as Hubble for 
being the discoverer of the expansion of the universe [J, Q| ■ 



In 1930, Eddington, who was working with McVittic 
on possible dynamical solutions of the cosmological Ein- 
stein equations (following a discussion that he had with 
de Sitter during a meeting of the Royal Astronomical 
Society [2l|), was alerted by Lemaitre about his own re- 
sults. He publicized them in a paper [22j in which he 
showed that the Einstein static universe is unstable (a 
result that was implicit in the paper of Lemaitre). He 
also emphasized a non-singular solution of Lemaitre in 
which the universe starts from the static solution of Ein- 
stein for t — > — oo, then expands (this model is often 
called the Eddington-Lemaitre model). In parallel, de 
Sitter advertised Lemaitre's discovery of an expanding 
universe in several publications [2^, |24[. Together with 
the supporting observations of Hubbl e 117 1, a nd indepen- 
dent theoretical work by Robertson 19l 20j|. this really 
meant the demise of the static universe. 

In 1931, Einstein was forced to admit the theoretical 
and observational reality of the expansion of the uni- 
verse. He took an extre me p osition and renounced to 
the cosmological constant [25[ that he considered being 
the "biggest blunder" of his life (in contrast, Lemaitre [J] 
always defended the cosmological constant as one of the 
most fundamental ingredients of modern cosmology). In 
1932, Einstein published with de Sitter a one-page paper 
[2r| in which they consider an expanding universe in a 
Euclidean space (zero curvature k = 0) without pressure 
(p = 0) and without cosmological constant (A = 0). This 
is the celebrated Einstein-de Sitter (EdS) solution that 
describes the matter era. The main motivation of these 
authors was to show that an expansion of the universe is 
possible even without the cosmological constant. 

In the meantime, Lemaitre had gone one step further. 
If the universe is expanding, it may have started from a 
singularity in the past. In 1931, he proposed a model (27j 
in which the universe begins its expansion from a singular 
initial stat43, or more precisely, from a near-singularity 
that he called the "primeval atom" [29pl. This expan- 
sion is followed by a phase of stagnation during which 
the universe remains close to the static Einstein solution. 
In Lemaitre's model, it is during this phase of stagnation 
that galaxies form. At later times, the universe enters in 
the de Sitter phase of exponential expansion. A new con- 
troversy started with Einstein who did not believe in the 
singular model of Lemaitre that implies that the universe 
has a "beginning" . This controversy remained until the 
end of their life. On the other hand, in 1960, Fred Hoyle, 



3 Eddington 1221 also found singular solutions to the cosmological 
equations but rejected them as being unphysical. In 1931, he 
wrote: "Philosophically, the notion of a beginning of the present 
order of Nature is repugnant to me" [28| . 

4 Lemaitre |30| also considered a model of universe in perpetual 
oscillation, undergoing cycles of expansion and contraction. He 
poetically called it "phoenix-universe". However, he finally re- 
jected this model at regret because he found that the duration 
of a cycle is too short. Cyclic universes were also discussed by 
Fricdmann 0, [j| and Einstein [3^] . 
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who promoted a "steady state theory" [3l| in which the 
universe has always been identical to itself and in which 
matter is created spontaneously, and continuously, made 
fun of Lemaitre during a meeting in Pasadena by intro- 
ducing him as the "Big Bang man" This is how the 
term "Big Bang" was introduced in cosmology and be- 
came popular. 

Up to the mid-1960's, it was not clear whether the 
early universe had been hot or cold. The discovery by 
Penzias and Wilson [32| of the 2.7 K cosmic microwave 
background radiation (CMB), whose existence had been 
predicted by the hot universe theory HH-HU , immediately 
led widespread acceptance to the Big Bang theory. Ac- 
cording to that theory, the early universe was filled with 
an ultrarelativistic gas of photons, electrons, positrons, 
quarks, antiquarks, etc. This is described by an equation 
of state p = pc 2 /3 from which it is easily deduced that the 
scale factor evolves as a oc t 1 / 2 , the density as p oc £~ 2 , 
and the temperature as T oc t~ x l 2 . At t = 0, the scale 
factor vanishes while the energy density and the temper- 
ature become infinite. This is why the point t = is 
known as the point of the initial cosmological singularity 
(Big Bang). As the universe expands, the density and 
the temperature decrease. When the universe is suffi- 
ciently cooled, particles and antiparticles annihilate each 
other, the photon-gas energy density falls off relatively 
rapidly and the main contribution of the matter density 
starts to come from the small excess of baryons over an- 
tibaryons, as well as from other fields and particles which 
now comprise the dark matter. 

However, the hot universe theory suffers from a se- 
ries of difficulties. The first difficulty is the singular- 
ity problem. The energy density goes to +oo at t = 
and the solution cannot be formally continued for t < 0. 
A natural question is: "What was before t = 0?" . Of 
course, one can argue that space and time started with 
the Big Bang, so that this question is meaningless. How- 
ever, this answer is not quite satisfactory. On the other 
hand, close to the Big Bang, the density can be arbi- 
trarily large. This is in conflict with the laws of physics. 
When the density becomes high enough, quantum effects 
should be taken into account in the problem^. As a re- 
sult, the Planck density p P = c 5 /G 2 h = 5.16 10" g/m 3 , 
the Planck length l P = (GTi/c 3 ) 1 / 2 = 1.62HT 35 m, the 
Planck mass M P = {fic/G) 1 ' 2 = 2.171CT 5 g, the Planck 
temperature T P = M P c 2 /k B = 1.42 1 32 K, and the 
Planck time t P = (hG/c b ) 1 / 2 = 5.3910~ 44 s should play 
a fundamental role in the very early universe. 

The standard Big Bang theory suffers from other prob- 
lems that are known as the flatness problem (or the fine- 
tuning problem), the horizon problem (or the causality 



a Remarkably, the question of the origin of space and time, and 
the importance of quantum mechanics in the early universe were 
discussed early by Lemaitre l36l : "If the world has begun with a 
single quantum, the notion of space and time would altogether 
fail to have any meaning at the beginning" . 



problem), and the monopole problem. These problems 
can be solved at one stroke in the framework of infla- 
tionary cosmology [13, [Hi- In that scenario, the universe 
was in an unstable vacuum-like state described by some 
hypothetical scalar field <p with its origin in quantum fluc- 
tuations. This leads to an equation of state p = —pc 2 , 
implying a constant energy density p, of the order of 
the Planck density p P , called the vacuum energy. As a 
result of this vacuum energy, the universe expands expo- 
nentially rapidly in a tiniest fraction of a second. This 
"inflation" can solve the above-mentioned problems. At 
the same time, the inflation theory predicts the scale in- 
variant spectrum of the primordial density fluctuations 
and explains why the universe is spatially flat (or be- 
haves as if it were spatially flat due to the smallness of 
k/a 2 ). This is therefore a very attractive theory. 

In this paper, we propose a model of the early universe 
based on a polytropic equation of state with index n = 
1. Before that, we first explain how we came to such a 
description. 

In previous works (39rj4l| . we explored the possibil- 
ity that the universe is filled with self-gravitating Bosc- 
Einstcin condensates (BEC) with short-range interac- 
tions (see a detailed list of references in (3(| and recent 
additional references in [42j). In the strong coupling limit 
(the so-called Thomas- Fermi approximation), the equa- 
tion of state of the BEC is p = (2ira s h 2 /m 3 )p 2 , where 
a s is the scattering length J43[. This is the equation of 
state of a polytrope with index n = 1 (TJ . Initially, 
this model was introduced in order to describe dark mat- 
ter halos. A nice feature of the BEC model is that it 
avoids the presence of cusps at the center of the halos 
because of the effective repulsion due to the Heisenberg 
principle or because of the self-interaction of the bosons 
which plays the same role as the Pauli exclusion principle 
for fermionf0. Since cusps are not observed in real dark 
matter halos, the BEC model (or the fermion model) is 
favored as compared to the standard cold dark matter 
(CDM) model that predicts cuspy density profiles. Of 
course, at large scales, these models become equivalent. 

The BEC model was then applied to cosmology in or- 
der to describe the evolution of the universe as a whole 
[45l-f50l] . In [13 , we investigated the evolution of a "BEC 
universe" with an equation of state of the form p = kp 2 c 2 . 
We considered the case of repulsive (k > 0) and attrac- 
tive (k < 0) self-interaction. More generally, we studied 
a "cosmic fluid" described by an equation of state of the 
form p = (ap + kp 2 )c 2 . When k = 0, we recover the 
standard linear equation of state p = ape 2 describing 
radiation (a = 1/3), pressureless matter (a = 0), stiff 
matter (a = 1), and vacuum energy (a = —1). There- 
fore, the generalized equation of state p = (ap + kp 2 )c 2 
can be viewed as the sum of a linear term describing 



6 If dark matter halos are made of self-gravitating fermions, the 
cusps are prevented by the Pauli exclusion principle; see |3S|| for 
discussions and references. 
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a "classical" universe and a quadratic term due to the 
BEC. At late times, when the density is low, the effect 
of the BEC is negligible and we recover the classical uni- 
verse. On the contrary, the effect of the BEC becomes 
important in the early universe where the density is high. 
Therefore, it is interesting to study how the BEC affects 
the early evolution of the universe. In our previous paper 
[47| , we assumed that BECs form after the radiation era, 
when the temperature has sufficiently decreased. This 
prevented us from extrapolating the solution before the 
radiation era (see the Remark in Sec. 8.2 of |47|). In the 
present paper, we explore the possibility that the equa- 
tion of state p = kp 2 c 2 may hold before the radiation era. 
Therefore, we assume that the early universe is described 
by an equation of state of the form p = (ap + kp 2 )c 2 with 
a = 1/3. Actually, the study of this generalized equation 
of state is interesting even if its origin is not connected to 
BECs. Furthermore, for this class of models, the Fried- 
mann equations can be solved analytically, which is an 
additional source of motivation. First results were given 
in [H, |47} where it was shown that the repulsive and 
attractive models behave very differently. For repulsive 
self- interactions (k > 0), the universe starts at t = from 
a singularity at which the density is infinite but the ra- 
dius finite. For attractive self- interactions (k < 0), the 
universe always existed in the past. For t — > —00, its 
density tends to a constant value while its radius tends 
to zero exponentially rapidly. 

In the present paper, we complete the study of the 
generalized equation of state p = (ap + kp 2 )c 2 . First, 
we analyze its thermodynamical properties and study 
the temporal evolution of the temperature. Secondly, 
for k < 0, we show that the exponential growth of the 
scale factor can account for a phase of inflation with con- 
stant density that we identify with the Planck density 
pp. Finally, taking a = 1/3, we study the transition 
between a pre-radiation era and the radiation era. We 
mention the analogy with a second order phase transi- 
tion where the control parameter is the time t, the order 
parameter is the scale factor a(t), and the Planck con- 
stant h plays the role of finite size effects (the standard 
Big Bang theory is recovered for h = 0). For k < 0, the 
radius of the universe passes from a, = lp = 1.62 10~ 35 m 
at t = U = to ai = 2.61 10~ 6 m at t = h = 1.25 10~ 42 s 
before entering in the radiation era and increasing al- 
gebraically as t x l 2 . During the inflationary phase, the 
density remains approximately constant with the Planck 
value (p ~ pp = 5.16 10"g/m 3 ) before decreasing as 
t~ 2 . The temperature passes from T 4 = 5.5410~ 173 K at 
t = U = to Ti = 3.93 10 31 K at t = h before decreasing 
algebraically as ~ t -1 / 2 . It achieves its maximum value 
T e = 7.4010 31 K at t = t e = 1.271(r 42 s. For k < 0, 
the initial radius of the universe a(0) = a\ = 2.61 10 -6 m 
is already "large" so there is no inflation. The veloc- 
ity of sound is less than the speed of light only for 
t>U = 8.32 10~ 45 s, corresponding to a; = 3.90 lCT 6 m, 
Pi = 1.2910"g/m 3 , and T t = 1.64 1 32 K. This marks 
the beginning of the physical universe in this model. 



The paper is organized as follows. In Sec. [HI we re- 
call the basic equations of cosmology and discuss famous 
cosmological models. In Sees. IIIII and IIV1 we study the 
generalized equation of state p = (ap + kp 1+1 / n )c 2 for 
any value of the parameters a, k and n > 0. We discuss 
two classes of solutions depending whether k is positive or 
negative. In Sees. IVland lVll we consider a specific model 
of physical interest corresponding to a = 1/3, n = 1, and 
k = ±4/(3pp). For k < (negative polytropic pres- 
sure), we get a model of inflationary universe without 
singularity. It describes in a unified manner the transi- 
tion from a pre-radiation era to the radiation era. For 
k > (positive polytropic pressure), we get a model of 
non-inflationary universe with a new type of initial sin- 
gularity. In Sec. IVII| we develop an analogy with second 
order phase transitions where the Planck constant plays 
the role of finite size effects. 



II. BASIC EQUATIONS OF COSMOLOGY 

In this section, we recapitulate the basic equations of 
cosmology [5l[ and discuss some famous models of uni- 
verse that were mentioned in the Introduction. This will 
prepare the ground for the next sections in which we dis- 
cuss generalized cosmological models. Some readers may 
directly go to Sec. IIIII 

A. The Einstein equations 

In a space with uniform curvature, the line element 
is given by the Friedmann-Lemaitre-Roberston- Walker 
(FLRW) metric 

ds 2 = c 2 dt 2 - a(t) 2 i^r^ + r 2 (d9 2 + sin 2 6 d<t> 2 )\ , 
L 1 — kr 2 J 

(1) 

where a(t) represents the radius of curvature of the 3- 
dimensional space, or the scale factor. By an abuse of 
language, we shall call it the "radius of the universe". 
On the other hand, k determines the curvature of space. 
The universe can be closed (k > 0), flat (k = 0), or open 
(k < 0). 

If the universe is isotropic and homogeneous at all 
points in conformity with the line element ([T]), and con- 
tains a uniform perfect fluid of energy density e(t) = 
p(t)c 2 and isotropic pressure p(t), the energy-momentum 
tensor is 

T ° = pc 2 , Tl =T 2 = Tl = -p. (2) 
The Einstein equations 

^ (3) 

relate the geometrical structure of the spacetime (gij) to 
the material content of the universe (T^). For the sake 
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of generality, we have included the cosmological constant 
A. Given (fTJ) and ©, these equations reduce to 



8irGp + A = 3- 



kc 2 



iirG lad - 
— p - A = 



(4) 



(5) 



where dots denote differentiation with respect to time. 
These are the well-known cosmological equations first de- 
rived by Friedmann [|| Q . 



the very early universe in which the pressure of radiation 
dominates over the dark energy (cosmological constant). 
On the other hand, we consider a flat universe (k = 0) in 
agreement with the observations of the cosmic microwave 
background (CMB) [52[. Then, the Friedmann equations 
reduce to 



dp 
~dt 

a 
a 



0, 



4ttG 



3p 



(12) 



(13) 



B. The Friedmann equations 

The Friedmann equations are usually written in the 
form 



d P , q d ( L P I 



a 
a 



H 2 = 



AttG 



8irG 



A 

3' 



cc 
7T 



A 

3' 



(6) 



(7) 



(8) 



where we have introduced the Hubble parameter H = 
a I a. Among these three equations, only two are inde- 
pendent. The first equation, which can be viewed as an 
"equation of continuity" , can be directly derived from the 
conservation of the energy momentum tensor d{T %3 = 
which results from the Bianchi identities. For a given 
barotropic equation of state p = p(p), it determines the 
relation between the density and the scale factor. Then, 
the temporal evolution of the scale factor is given by Eq. 



Equivalent expressions of the "equation of continuity" 



arc 



da 



(pa 3 ) 



nP 2 

& 



dt dt 1 Kl 



pc 2 )} ■ 



(9) 



(10) 



(11) 



Introducing the volume V = (4/3)7ra 3 and the energy 
E = pc 2 V : Eq. (fTQ|) becomes dE = —pdV. This can be 
viewed as the first principle of thermodynamics for an 
adiabatic evolution of the universe dS = 0. This thermo- 
dynamical interpretation of the Friedmann equation (|6]) 
was first given by Lemaitre [TH, HH . 

In most parts of this paper (except in Sec. Ill D|) . we 
shall take A = 0. This is because we are interested in 



# = I - I = 

a , 



8ttG 



(14) 



The Friedmann equations can also be derived in the 
context of Newtonian cosmology. In this framework, 
Milne [H| and McCrea & Milne [H] first obtained the 
Friedmann equations without pressure (i.e. for c — >■ 
+oo). Later, McCrea [55| introduced relativistic effects 
in Newtonian cosmology and obtained the complete set of 
Friedmann equations with pressure. Another derivation 
was given by Harrisson [56 1 who obtained a relativistic 
generalization of the Euler-Poisson system in the context 
of Newtonian cosmology. Finally, the continuity equation 
was corrected in a recent paper by Lima et al. |57| . Actu- 
ally, the derivation of the Friedmann equations with rela- 
tivistic effects (pressure) from Newtonian cosmology was 
first exposed by Lemaitre [29| as early as 1931 (see also 
[jj|). However, for Lemaitre, the Newtonian derivation 
of the cosmological equations is just a pedagogical means 
and does not reflect the true (relativistic) nature of the 
problem. In modern textbooks of astrophysics [H, [59| , 
Newtonian cosmology with pressure is presented as a rig- 
orously valid description of the universe in a region small 
compared to the Hubble length c/H. 

Remark: Considering Eq. (J7JI, we see that the cosmo- 
logical constant corresponds to a force F = Aa/3 in the 
Newtonian interpretation. A positive cosmological con- 
stant (de Sitter) acts as a force of cosmic repulsion and a 
negative cosmological constant (anti de Sitter) acts as a 
force of cosmic attraction whose effect is similar to that of 
a harmonic oscillator. On the other hand, a positive pres- 
sure has an attractive character while a negative pressure 
has a repulsive character. Gravity is attractive (the uni- 
verse decelerates) when p + 3p/c 2 > and repulsive (the 
universe accelerates) when p + 3p/c 2 < 0. 



C. The thermodynamical equation 



From the first principle of thermodynamics 

dS(V, T) = i [d(p(T)c 2 V) + P (T)dV] , (15) 



one can derive the thermodynamical equation [51j : 

-Xpc 2 +p). (16) 



dp 1, 2 



dT T v 



6 



For a given barotropic equation of state p = p(p), this 
equation can be used to obtain the relation T = T(p) 
between the temperature and the density. Combining 
Eq. (HU) with Eq. (pi]), we get 



dS_ 
~dt 



0. 



where 



S=^(p + pc 2 ). 



(17) 



(18) 



is the entropy of the universe in a volume a 3 . This con- 
firms that the Friedmann equations (H])-© imply the con- 
servation of the entropy [5l|. 



a,ds = c-\/3/A is the radius of the original (static) de Sit- 
ter model. The solution (|21|) describes an empty universe 
(p = 0) that expands because of the repulsive nature of 
the cosmological constant. Any model with A > goes 
to a de Sitter model for a — > +oo. 

The Lemaitre model of 1927 cor resp onds to k = +1, 
A > 0, p = and p > [H The equation 

of continuity ^ implies pa 3 oc 1 which can be inter- 
preted as the conservation of mass-energy for a pressure- 
less universe. We can write this relation in the form 
p = aM e I '(2ir 2 'a 3 ) , where Me is Einstein's mass and a 
is a dimensionless constant. Introducing the dimension- 
less scale factor R = a/aE and the dimensionless time 
we can rewrite the Friedmann equation (|8|) in 

the form 



D. Famous models of the universe 

In 1917, Einstein Q constructed a static model of the 
universe. He assumed that the universe if filled with a 
pressureless gas {p = 0) and introduced a cosmological 
constant A in the equations of general relativitjQ- Setting 
d = a = in Eqs. it is necessary that A > and 

k = +1. Then, 



a E 



c 



PE 



A 

4^G 



47rGa| ' 



(19) 



The static Einstein universe is finite (though unbounded) 
with a positive curvature and a density fixed by A and 
G. The volume of this universe is Ve = 2ir 2 a E - Its mass 
Me = PeVe is therefore 



M f 



7TC 



2GVA 



TTC 2 aE 

2G 



2tt 



2 (— ) 



3/2 



1 



(20) 



The same year, de Sitter found another cosmolog- 
ical solution of the Einstein field equations corresponding 
to an empty universe (p = p = 0) with a cosmologi- 
cal constant A > 0. In his original article, he aimed 
at constructing a static solution. He obtained a radius 
ads = Cy/3/A. However, it was understood later by 
Lemaitre [14j that, by using a proper change of variables 
in the metric, the de Sitter solution actually describes a 
dynamical universe whose radius increases exponentially 
rapidly with time. Setting p = p = and k = in Eqs. 
©-©, we get 



a(t) = a(0)e^ 



(21) 



In this model, the Hubble parameter is constant: H = 
a I a = ^/A/3. It can be written H = c/aas where 



7 Actually, there exist another static solution of the Einstein 
equations without cosmological constant corresponding to p = 
-(l/3)pc 2 , k = +1 and p = 3c 2 /(8vrGa 2 ). 



R\ 2a 11 
R ) ~3R 3 ~R 2 + 3' 



(22) 



For a = 1, this equation admits the static Einstein solu- 
tion R = 1 (i.e. a = aE)- However, as shown by Lemaitre 
(l5l [HI , and emphasized by Eddington [22j , this solution 
is unstable and, for sufficiently large times, the radius of 
the universe expands. This solution, according to which 
the radius of the universe is equal to the Einstein radius 
for t — > — oo and increases for large times (behaving ul- 
timately like the de Sitter solution (f2"Tj) for r — > +oo) is 
called the Eddington-Lcmaitre model (see Fig. [1]). This 
model does not have an initial singularity. 

For a = 1, the differential equation (|2"2"j) has the ana- 
lytical solution 



r=(R-l] 



R 



R 3 -3R + 2 



2\/3 Argsinh 




+ ln 



Vr- i 
i+i 



In 



^3(2 + R) + 2- Vfl 
y/3(2 + R) + 2 + Vr 



where r is defined up to an additional constant. 
t — > — oo, we obtain 



R-l 



-Ci 



(Einstein) 



(23) 



For 



(24) 



with Ci = 2y/3 Argsinh(l/V2) - ln(6) 
t — > +oo, we get 



0.489279. For 



R 



1 



(r+C 2 )/V3 



(de Sitter), 



(25) 



with C 2 = ln[(V3 + 1)/(V3 - 1)] ~ 1.31696. To the best 
of our knowledge, these expressions have not been given 
previously. 

In 1931, Lemaitre |27| considered a model in which a 
is slightly larger than one (i.e. it contains more mat- 
ter than the static Einstein model). This model starts 
from a singularity at t = 0. The radius first increases 
as R ~ (3a/2) 1 / 3 r 2 / 3 (this corresponds to the solution 
(|3"3"|) found later by Einstein & de Sitter [26[), then slows 
down and reaches a minimum rate at R ~ a 1 / 3 , after 
which it re-expands rapidly, ultimately approaching the 



7 



Lemaitre models 



4- 



de Sitter 



2 - 



Einstein 



Stagnation period 




with — 1 < a < 1. The pressure is positive for a > 
and negative for a < 0. For this equation of state, the 
Fricdmann equations (fT2|) leads to the relation 



pa 



3(l+a) 



CX 1. 



(29) 



In the following, we assume (the case a = — 1 is 

treated in Sec. IIIE4j) . Substituting Eq. (J29]) in Eq. (fT4|) 
and solving the resulting equation for a(i), we find that 
the radius of the universe increases with time as 



aKt 2/[3(l + a)]_ 



(30) 



FIG. 1. Evolution of the scale factor in the Eddington- 
Lemaltre universe (dashed-line) and in the Lemaitre universe 
(full line). We have taken a = 1.00001. 



We have considered expanding solutions and defined the 
origin of time t = such that a(0) = 0. The Hubble 
parameter and the density decrease with time as 



a 

H = - = 



a 3(1 + a)t ' 



6ttG(1 + a)H 2 ' 



(31) 



de Sitter solution (f2"T|) for r — > +oo. The main feature 
of this model is the existence of a "stagnation period"@ 
during which the radius remains close to R ~ a 1 / 3 . Dur- 
ing the stagnation period, the differential equation (f2"2")l 
takes the approximate form 



In all these models, the universe is singular at t = 0. 
The radius vanishes while the density and the Hubble 
parameter are infinite (Big Bang singularity). 

For the equation of state (j2"8")l . the thermodynamical 
equation (fl"6]l can be integrated into 



whose solution is 1511: 



R 



!/ 3 fl + (l-a- 2 /Y /2 smh(r-T m ) 



(26) 



(27) 



rocp «/(a+l) __L 
- 3a 



(32) 



where r m is the time at which R reaches its minimum. If 
a is close to unity, then a will remain close to the Einstein 

-2/3 



The temperature decreases with time when a > and 
increases with time when a < 0. 

According to Eq. (fl~3"|) . the universe is decelerating for 
a > —1/3 and accelerating for a < —1/3. For a = —1/3, 
the scale factor increases linearly with time (a cut). 
Finally, the square of the velocity of sound is c? s = 



static value a E for a long time At — | ln(l — a ' )|. p'(p) = ac 2 . The velocity of sound is real for a > (it is 



For a = 1, we recover the Eddington-Lemaitre model. 
However, the historical interest of the Lemaitre model 
[l?! of 1931 is that it predicts for the first time that the 
universe may have emerged from a singularity (or near- 
singularity) at t = in which the radius is equal to zero 
and the density is infinite. This is what Lemaitre called 
the "primeval atom" [2Sj . and which became later the 
Big Bang theory. 



E. Linear equation of state 

In early works of cosmology, it was assumed that the 
universe is pressureless (p = 0). In later works, pressure 
effects have been taken into account. Many workers have 
considered a linear equation of state 



less than the speed of light for a < 1) and imaginary for 
a < 0. It vanishes for a = 0. 



1. Matter ( Einstein- de Sitter universe) 

A pressureless universe (p = 0) made of baryonic mat- 
ter or dark matter (dust) corresponds to a = 0. This is 
the so-called Einstein-de Sitter (EdS) universe. In this 
model, 



2/3 



p = ape , 



(28) 



p cx 1/a 3 , oat 



H=- = J_ 

3*' 9 6nGt 2 ' 



2. Radiation 



(33) 



(34) 



The equation of state of a radiation (photon gas) is 

Lemaitre thought that galaxies formed during this stagnation V = pc 2 /3, corresponding to a = 1/3. This yields 



period, and that this "instability" triggered the expansion of the 
universe. 



pcxl/a 4 , acxt 1/2 , (35) 
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T 1/4 1 1 
T oc p 1 oc - oc -— . 

a t 1 / 2 



(36) 



(37) 



The relation between the energy density and the temper- 
ature corresponds to the Stcfan-Boltzmann law 



pc 2 = oT\ 



15c 3 h 3 ' 



(38) 



Introducing the Planck scales, we can write the fore- 
going equations in the form 



3pp ( tp 



T = Tp 



T = Tp 



32tt V t 



\ 1/4 (_P_^ 



45 \ 1/4 (tp^ 112 



32tt s 



(39) 



(40) 



(41) 



3. Stiff equation of state 

A stiff equation of state p = pc 2 , for which the velocity 
of sound c s = (dp/dp) 1 ' 2 is equal to the speed of light c, 
corresponds to a = 1. In that case, 



interpreted p\ as the vacuum energy. The origin of the 
vacuum energy was later discussed by Sakharov [6(3| and 
Zeldovich [6l| in relation to quantum field theory. How- 
ever, there is a fundamental difficulty in interpreting p\ 
as the vacuum energy because quantum field theory pre- 
dicts that the vacuum energy density should be of the 
order of pp which is 122 orders of magnitude larger than 
Pa- This is known as the cosmological constant problem 
[62]. Therefore, p\ may not be related to the vacuum 
energy (hence to quantum fluctuations) but could be a 
fundamental property of general relativity (see discussion 
in Paper II). In the modern terminology, the effect of the 
cosmological constant, or the effect of an exotic fluid, on 
the accelerated expansion of the universe is called dark 
energy [63j . 

The equation of state p = —pc 2 also appeared in the 
paper of McCrea 1551 in his discussion of the steady state 
theory of Hoyle [31|. In this theory, in which there is a 
continuous creation of matter, the density and the Hub- 
ble parameter are constant, leading to Eq. (|4l)|) . Accord- 
ing to Hoyle [3l[ , a continuous creation of matter can be 
achieved by a proper modification of the Einstein equa- 
tions. McCrea [55| noted that the same result could be 
obtained with an equation of state p = —pc 2 without in- 
voking creation of matter and without having to modify 
the Einstein equations. 

The equation of state p = —pc 2 als o ap peared in in- 
flationary models of the early universe [38| . It leads to a 
constant density which is identified with the Planck den- 
sity pp. In that case, this equation of state may well be 
related to the vacuum energy density. 



p oc 1/a 6 , a oc t 1 ^ 3 , 



(42) 



F. Necessity of the inflation 



H 



1 



3t' 



2AnGt 2 '' 



ml/2,11 

1 oc p 1 oc — n oc — . 
1 a 3 t 



(43) 



(44) 



4- Vacuum energy and dark energy 

For the equation of state p — —pc 2 , corresponding to 
a = — 1, the density and the Hubble parameter are con- 
stant 



p oc 1, 



H 



lirGp 



(45) 



The previous models, based on a linear equation of 
state, lead to a finite time singularity at t = (except 
the model with a — —1). This poses several problems. 
We briefly recall what the problems are and show in the 
next sections how a modification of the equation of state 
can resolve them. 

We call Ho the present value of the Hubble constant 
(the subscript refers to present-day values). The Hub- 
ble time is I/Hq. It gives an estimate of the age of the 
universe. The characteristic "size" of the universe (dis- 
tance of the cosmological horizon) is the Hubble length 
a Q = c/H . Using H = 2.27 1CT 18 s _1 , we find that 
1/H = 4.41 10 17 s and a = 1.32 10 26 m. For a flat space 
(k = 0), the density of the universe po can be obtained 
from the Friedmann equation 



and the scale factor increases exponentially rapidly as 



a(t) = o(0) 



Ht 



(46) 



Lemaitre |58( was the first to understand that the ef- 
fect of the cosmological constant A is equivalent to that 
of a fluid with density p = p\ = A/8irG [see Eq. (0) or 
HJ] described by an equation of state p = —pc 2 . He 



po 



3H 2 _ _3_ 
8ttG ~ 8tt 



(H Q tp) 2 p P . 



(47) 



Numerically, po = 9.20 1 24 g/m 3 . On the other hand, 
the temperature measured from the cosmic microwave 
background (CMB) is T = 2.73 K. 

The early universe is dominated by radiation. In the 
standard model, the radiation phase leads to a finite time 
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singularity at t = (Big Bang). It is expected that 
quantum mechanics will come into play at the Planck 
time tp when the universe is very "small". At that time, 
the density and the temperature are (see Sec. Ill E 2|) : 



p(tp) = — pp = 1.5410 98 g/cm 3 , 



and 



T{t P ) = 



( 45 



V32tt j 



1/4 



T P = 6.5410 J1 K. 



(48) 



(49) 



The particle horizon at the Planck time is of the order 
of the Planck length l P = ct P = 1.621CT 35 m. Since 
T{t)a{t) is invariant during the evolution of the universe, 
we have 



45 
32V 



1/4 



T P a(t P ) = T a . 



(50) 



If we take a(tp) = lp, we find that the present size 
of the universe is of order 3.8710 _4 m which is obvi- 
ously too small. The same argument can be turned the 
other way round. Using the present size of the universe 
ao = 1.32 1 26 m, we find that its size at the Planck time 
is a(tp) = 5.5210~ 6 m which is much larger than the 
Planck length (p by a factor ~ 10 29 . 

It is instructive to recover this result in a slightly 
different, but equivalent, manner. Since p ra d{t)o,(t) i s 
conserved during the evolution of the universe, we have 
p{tp)a{t P ) 4 = Prad,oa^. Writing p rad ,o = ^rad,oPo, and 
using Eq. (|47|) . we obtain 



a(tp) 
ao 



= {^rad,o) 1/4 VHotp = 4.75 10" 



Since 



no 



H Q t P = 1.2210 



-61 



the foregoing equation can be rewritten 



a(tj 



h 



] - = (4O. a d, ) 1/4 -ri 



VHoTp~ 



= 3.8810 



oil 



(51) 



(52) 



(53) 



Since H tp = 1.22 10 -61 (meaning that the Hubble time 
1/Hq is tremendously larger than the Planck time tp), 
we see that 



lp <C a(t P ) ^ a . 



(54) 



The fact that a{tp) ^s> lp makes problem [38[. One must 
find a mechanism to bring the "size" of the universe from 
l P = 1.62 10~ 35 m to a(t P ) = 6.28 10~ 6 cm in a very short 
lapse of time. This can be achieved with the inflationary 
scenario. The usual inflationary scenario is based on the 
equation of state p = —pc 2 (vacuum energy). As we 
have seen in Sec. Ill E 41 this implies that p and H are 
constant so that the scale factor increases exponentially 



rapidly with time, as e Ht with H = (8nGp/3) 1/2 . If 
we assume that p ~ pp, we obtain H = (87r/3) 1 ' 2 tp 1 
so that the exponential time scale is of the order of the 
Planck time. This solves the above-mentioned problem 
[38l |. In the following sections, we show that a poly tropic 
equation of state with n > produces similar results and 
leads to simple cosmological models unifying the phases 
of inflation and radiation. 



III. GENERALIZED EQUATION OF STATE 
WITH POSITIVE INDEX 



We consider a generalized equation of state of the form 

p={ap + kp 1+1/n )c 2 . (55) 

This equation of state was introduced in our previous 
paper [47[ for the specific index n = 1. This is the sum 
of a standard linear equation of state p = ape 2 and a 
polytropic equation of state p = kp^c 2 , where k is the 
polytropic constant and 7 = 1 + 1/n is the polytropic in- 
dex. Concerning the linear equation of state, we assume 
— 1 < a < 1 (the case a = — 1 is treated specifically in 
Appendix [XJ ■ This equation of state describes radiation 
(a = 1/3), pressureless matter (a = 0), and vacuum en- 
ergy (a = —1). Concerning the polytropic equation of 
state, we remain very general. This polytropic equation 
of state may correspond to self-gravitating BECs with re- 
pulsive (k > 0) or attractive (k < 0) self-interaction (the 
standard BEC model corresponds to n = 1) [471 ] . but it 
may have another origin. 

For an equation of state of the form ([55]) with posi- 
tive index n > 0, the polytropic component dominates 
the linear component when the density is high. These 
models, studied in the present paper, describe the early 
universe. Conversely, when n < 0, the polytropic compo- 
nent dominates the linear component when the density is 
low. These models, studied in Paper II, describe the late 
universe. As will be discussed in detail in Paper II, these 
two situations are strikingly similar. They reveal a form 
of "symmetry" between the early and the late universe. 
In Papers I and II, we assume that a + 1 + fcp 1 /" > 0. 
This corresponds to the "normal" case where the den- 
sity decreases with the radius (see below). The opposite 
case a + 1 + fcp 1 /" < 0, leading to a "phantom universe" 
where the density increases with the radius, is considered 
in Paper III. 

The purpose of our series of papers is to make an ex- 
haustive study of all the possible cases (a,n,fc), even if 
some of them seem to be unphysical or in conflict with 
the known properties of our universe. Actually, by re- 
jecting the cases leading to past or future singularities, 
we will be led naturally to the "good" model (see Paper 
II). 
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A. The density 



B. The temperature 



For the equation of state (|55[) . the Friedmann equation 
(Ti"2"j) becomes 

^ + 3-p(l + a + kp^ n ) =0. (56) 
at a 

Assuming a + 1 + hp 1 / 71 > 0, this equation can be inte- 
grated into 



(57) 



[(a/a*) 3 ( 1+Q )/" + 1] 



where p* = [(a + l)/|fc|]" and a* is a constant of inte- 
gration. The upper sign corresponds to fc > (repulsive 
self-interaction) and the lower sign corresponds to fc < 
(attractive self- interaction) . We may draw a parallel be- 
tween these two types of distributions and the statistics 
of fermions (+) and bosons (— ). 



1.5 



0.5 





1 i | 

k>0 


k<0 








i 


• i 



2 

a/a. 



FIG. 2. Density as a function of the scale factor for k > 
and k < 0. We have taken n = 1 and a = 1/3. 



For k > 0, the density is defined only for a > a*. When 
a — > a*, 



p* 



3(1 + a) 



1 



(a/a* - 1)" 



(58) 



When a — > +oo, p/p* ~ (a*/a) 3 ( 1+Q ) — > corresponding 
to the linear equation of state. In the same limits, p —> 
+oo and p — > 0, respectively. 

For k < 0, the density is defined for all a. When a -+ 0, 
the density tends to a finite value /?* When a — > +oo, p ~ 
/9*(a»/a) 3 ' 1+Q ' 1 — > corresponding to the linear equation 
of state. In the same limits, p — > —p*c 2 and p — > 0, 
respectively. 

Some curves giving the evolution of the density p as 
a function of the scale factor a are plotted in Fig. [2] for 
k > and fc < 0. 



For the equation of state (|55j) . the thermodynamical 
equation ([To]) becomes 

(" + fc7P 1/ ")^ = ^(a + l + fcp 1/ >. (59) 
This equation can be integrated into 



T = T» 



1 ± (P/P*) lfn 



(a+n+l)/(a+l) 



a/(a+l) 



(60) 

where T* is a constant of integration. This relation be- 
tween the temperature and the density can be viewed as 
a generalized Stcfan-Boltzmann law. Combined with Eq. 
(|57|) . we obtain 



T = T* 



(a/a,) 



3(a+n+l)/n 



[(a/a,) 3 (i+«)/" T l 



n+l 



(61) 



Let us consider some asymptotic limits. 
When a ^> a*, Eq. ([6T|) reduces to 

T/T, 



(a/a*) 3 "' 



(62) 



returning Eq. ([32]) for a linear equation of state. As we 
have already indicated, in this regime, the temperature 
decreases to zero when a > and increases to +00 when 
a < 0. When a = 0, it tends to a constant value T*. 
When fc > and a — > a*, 



T/T* 



n+l 



1 



(a/a* 



.3(1 + a) 
When k < and a -+ 0, 

T/T* ~ (a/a*) 3 ( Q+,i+1 )/' 1 



0. 



(63) 



(64) 



The extremum of temperature (when it exists) is lo- 
cated at 



Pe 
P* 



Te 
I 



+ 



a e ( a + n + 1 

— = + 

a* V an 



(a+l)(n + l). 

n/[3(a+l)] 



(65) 
(66) 



(a + l)(n+l) 



n + a + 1 



(l + a)(n+l) 



(67) 

Some curves giving the evolution of the temperature T 
as a function of the scale factor a are plotted in Figs. [3] 
and [4] for fc > and fc < 0. 

Finally, the entropy (fT8|) is given by 



(68) 



S = (a + l)-fp*c 2 , 
and we explicitly check that it is a constant. 
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FIG. 3. Temperature as a function of the scale factor for 
k > 0. We have taken n = 1, a = —1/3, a — and a = 1/3. 




a/a, 



FIG. 4. Temperature as a function of the scale factor for 
k < 0. We have taken n = 1, a = —1/3, a — and a = 1/3. 



FIG. 5. Evolution of w, q and (c s /c) 2 as a function of the 
scale factor for k > 0. We have taken n = 1 and a = 1/3. 
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- -(a+n+l)/n 
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FIG. 6. Evolution of w, q and (c s /c) 2 as a function of the 
scale factor for k < 0. We have taken n = 1 and a — 1/3. 



C. The equation of state parameter u>(i) 

We can rewrite the equation of state (|55j) as p 
w(t)pc 2 with 



w(t) = a± (a + 1) 



We define 



Pw _ ( a \ aw_ _ f 1_ 

p* \ a + 1 J ' a* V 01 



n/[3(a+l)] 



T 



(T«) 



an/ (a+l) 



(a 



i) 



n+l 



(69) 



(70) 



(71) 



corresponding to a possible point where the pressure van- 
ishes (w = 0). 

For fc > 0, to — >• +oo when a — > a, and u> — > a when 
a — > +oo. For a > 0, the pressure is always positive 
(u> > 0). For a < 0, the pressure is positive when a, < 
a < a w and negative when a > a w . 



For fc < 0, w — > — 1 when a — > and w — > a when 
a — > +oo. For a < 0, the pressure is always negative 
(w < 0). For a > 0, the pressure is negative when a < a w 
and positive when a > a w . 

Some curves giving the evolution of w as a function of 
the scale factor a are plotted in Figs. [5] and |B] for k > 
and fc < 0. 



D. The velocity of sound 

For the equation of state (|55p , the velocity of sound is 
given by 



P'(p) 



a±(a + l) 



n+l 



l/n 



(72) 



The velocity of sound vanishes at the point defined by 
Eqs. (|65|) -(|67 ]) where the temperature is extremum. At 
that point the pressure is extremum with value 



Pe 

p*c 2 



n+l 



+ 



(a + l)(n+l) 



(73) 
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The case c 2 < corresponds to an imaginary velocity of 
sound. We define 



Ps 

p* 



± 



(1 — a)n 



(l + a)(n + l) 



Os 

(I,: 



a + 2n + 1 



n(l — a) 



n/[3(l+a)] 



(74) 



(75) 



±- 



n(l 



a 



(a+l)(n + l) 



l + 2n 



a 



(1 + + 1) 



(76) 



corresponding to a possible point where the velocity of 
sound is equal to the speed of light. Different cases have 
to be considered. 

We first assume k > 0. When a — > a*, (c s /c) 2 — > +oo; 
when a — > +oo, (c s /c) 2 — > a. For a > 0, c 2 is always 
positive. For a < 0, c 2 is positive when a* < a < a e and 
negative when a > a e . The velocity of sound is larger 
than the speed of light when a* < a < a s and smaller 
when a> a s . 

We now assume k < 0. When a — > 0, (c s /c) 2 — > — (a + 
n + l)/n; when a — ^ +oo, (c s /c) 2 — > a. For a < 0, c 2 is 
always negative. For a > 0, c 2 is negative when a < a e 
and positive when a > a e . For a < 1, the velocity of 
sound is always smaller than the speed of light. 

Some curves giving the evolution of (c s /c) MS IX func- 
tion of the scale factor a are plotted in Figs. [5] and [6] for 
k > and k < 0. 



IV. EVOLUTION OF THE SCALE FACTOR 
A. The deceleration parameter 

The deceleration parameter is defined by 



, „ aa 
?*) = —• 



(77) 



The universe is decelerating when q > and accelerating 
when q < 0. Using the Fricdmann equations (|13[) and 
(fT4")) . we obtain for a flat universe 



<?(*) = 



1 + 3w(t) 



(78) 



For the equation of state using Eq. we get 



... l + 3a , 3, -t \ f P 

q{t) = — ± 2 {a+1) (i 



1/n 



(79) 



We define a critical density, scale factor, and temperature 
1 + 3a 



P* 



"3(1 + a) 



Oc = f 2 
a* V 1 + 3a 



n/[3(l+a)] 

(80) 



Tc _ (a+n+l)/(a+l) 



i)/( a +i) [=F(1 + 3a)] 



an/(a+l) 



[3(1 + a 



in+1 



(81) 



corresponding to a possible inflexion point (q = d = 0) 
in the curve a(t). Different cases have to be considered. 

We first assume k > 0. When a — > a*, g — > +oo; when 
a — > +oo, g — > (1 + 3a)/2. For a > —1/3, the universe is 
always decelerating (q > 0). For a < —1/3, the universe 
is decelerating when a* < a < a c and accelerating when 
a > a c . 

We now assume k < 0. When a — > 0, q — > —1; when 
a -+ +oo, q —> (1 + 3a)/2. For a < —1/3, the universe is 
always accelerating (q < 0). For a > —1/3, the universe 
is accelerating when a < a c and decelerating when a > 
a c . 

Some curves giving the evolution of q as a function of 
the scale factor a are plotted in Figs. [5] and |B] for k > 
and k < 0. 



B. The differential equation 

The temporal evolution of the scale factor a(t) is de- 
termined by the Friedmann equation (|14j) . Introducing 
the normalized radius R = a /a*, the density ([57| can be 
written 



P = 



[_R3(l+o)/n zp l]n ' 



(82) 



Substituting this expression in Eq. (|14j). we obtain the 
differential equation 



R = 



eKR 



[_R3(l+a)/n zp ljn/2' 



(83) 



where A' = (SttG/c^/S) 1 / 2 and e = ±1. In general, we 
shall select the sign e = +1 corresponding to an expand- 
ing universe (R > 0). The solution can be written as 



eKt 



/ 2 dR 



R 3(l+ a )/n- fl ( 84) 



R 



or, after a change of variables x — i? 3 ( 1+Q )/™ ; as 



3(a + l) 



eAt 



(x+1) 



n/2 



<7.r 



(85) 



The integral can be expressed in terms of hypergeomet- 
ric functions. Some simple analytical expressions can be 
obtained for specific values of n. Actually, we can have a 
good idea of the behavior of the solution of Eq. (|83|) by 
considering asymptotic limits (see below) . The complete 
solution is represented in the figures by solving Eq. 
numerically. 



C. The case k > 

The universe starts from a primordial singularity at 
t = (sec below for a revision of this statement) with a 
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finite radius R = 1 and an infinite density, pressure, and 
temperature. When t — > 0, 



i?~ 1 



P* 



3(1 + a 
3 



n/(n+2) 



2/(n+2) 



-2n/(n+2) 



,(86) 
(87) 



T 



3._, . n + 

- (1 + a) 

2 n 



-iff 



-2(n+l)/(n+2) 



Then, the universe expands indefinitely. When t — > +oo, 
the density decreases to zero and we recover the solution 
corresponding to a linear equation of state 

2/[3(l+a)] 



R 



P* 



-(l + a)Kt 



T 

T 



-{l + a)Kt 



-2q/(1+q) 



(89) 
(90) 
(91) 



For a > —1/3, the universe is always decelerating (i? < 
0). For a < —1/3, the universe is decelerating for R < 
R c = [-2/(l + 3a)] ,l /I 3 ( 1+Q) ] and accelerating for R > R c . 
The velocity of sound is less than the speed of light when 
R> R s = {(a + 2n + l)/[n(l - a)]}" /[3(1+a)1 . 
For n = 1, using the identity 

/dx 
Vx-1 — = 2Vx - 1 - 2 arctan (Vx- l) , (92) 

the general solution of Eq. 



is 

y/R3(i+<*) - 1 - arctan V R^ 1+a) - 1 = ^^"' Kt. 

(93) 

For a < —1/3, the time t c at which the universe starts 
accelerating is given by 



^±f?U. = t -3(1 + a) _ arctani /-3(1 + a) 



l + 3a 



1 + 3a 



The time f j at which the velocity of sound is equal to the 
speed of light is given by 



3(1 + a) 2(1 + a) 2(1 + a) , . 

>-KU = \ I \ ' - arctan \ I \ ' . (95) 



1-a 



1 - a 



This gives the time at which the physical universe begins. 
It may be meaningless to consider smaller times because 
the velocity of sound is greater than the speed of light. 
For n = 2, the general solution of Eq. (|53")) is 

fl*(i+")-|(a+l)lniZ= 3(1 + Q) i<rt + C. (96) 

Typical evolutions of the scale factor i?(t) correspond- 
ing to a < —1/3, a = 1/3, and a > 1/3 are represented 
in Fig. El 
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FIG. 7. Evolution of the scale factor in the case k > for 
different values of a (specifically a = —2/3, a — —1/3, and 
a — 1/3). We have taken n = 1. 



D. The case fc < 

The universe starts from t — > — oo with a vanishing 
radius R = 0, a vanishing temperature T = 0, and a 
finite density p = p*. When t — > — oo, 



R(x e 



Kt 



T 

T 



<-(a+n+l)Kt 



(97) 



This corresponds to an exponential expansion (early in- 
flation) due to the fact that the density is approximately 
constant. The universe expands indefinitely. When 
t — > +oo, the density decreases to zero and we recover 
the solution (|89|) -(f9"T |) corresponding to a linear equation 
of state. For a < —1/3, the universe is always accelerat- 
ing (R > 0). For a > —1/3, the universe is accelerating 
when R < R c = [2/(1 + 3a)] n ^ 3{ - 1+a ^ and decelerating 
when R > R c . In this model, the velocity of sound is 
always less than the speed of light. 
For n = 1, using the identity 



■ dx 



1— =2 Vx~+1 + Inf ill - ), (98) 
x Wx + 1 + 1/ 



(94) the general solution of Eq. (|8"3"|) is 



^/R3(i+») + l - In 



1 + Vi? 3 ( 1+a > + 1 



i?3(l+a)/2 

3(1 + a) 



Kt + C, (99) 



where C is a constant of integration. For a > —1/3, the 
time t c at which the universe starts decelerating is given 
by 



^±^Kt c + C= x m±^ 
2 V l + 3a 



h p + v^ ■ (100) 
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When t —> — oo, we obtain from Eq. (l99l) the asymptotic 
behavior 



Kt+D 



R(t) ~ e 

with D = 2(C + In 2 - 1)/[3(1 + a)]. 

For n = 2, the general solution of Eq. (|83p is 



(101) 



flf + -( a + l)lnR = 3(1 ^ a) Kt + C. (102) 

Typical evolutions of R(t) corresponding to a < —1/3, 
a = 1/3, and a > 1/3 are represented in Fig. [5J 
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FIG. 8. Evolution of the scale factor in the case k < for 
different values of a (specifically a = —2/3, a = —1/3, and 
a = 1/3). We have taken n = 1. 



A MODEL OF NON-SINGULAR 
INFLATIONARY UNIVERSE 



The generalized equation of state (55") with n > and 
fc < can be used to describe the transition between a 
phase of early inflation where the density of the universe 
is constant (p = p*) and a phase where the universe has 
a linear equation of state (j> = ape 2 ). Therefore, if we 
take a = 1/3, it describes the transition between the 
pre-radiation era and the radiation era. This provides 
a model of non-singular inflationary universe. For sim- 
plicity, we shall take n = 1 (in the BEC model, this 
corresponds to the standard index ji^). 



A. The basic equations 

For a = 1/3, n = 1, and k < 0, Eqs. (55), (57), (TO), 
and (6T) become 



p = (p/3 + kp 2 )c 2 , 



(103) 



(104) 



t = r, (l - /Vp*) 7/4 (p/ P *) 1/4 

(a/ ai ) 7 



T = T* 



[(«/«i) 4 + ir 



(105) 



(106) 



where p* = 4/(3|fc|), and we have noted a\ for a*. 

When oCai, the density has a constant value p ~ p*. 
This gives rise to a phase of early inflation. The value 
p* defines a fundamental upper bound p m ax for the den- 
sity. It is natural to identify it with the Planck density 
Pp = 5.16 10 99 g/m 3 (see Sec. Ill F[) . This is how quantum 
mechanics is taken into account in our model. Since the 
Planck density is obtained from a pure dimensional anal- 
ysis, we should write p* = upp, where k is a constant of 
order unity (actually, the Planck constant should be de- 
fined so as to correspond exactly to p m ax)- For simplicity, 
we shall take 



Pp- 



(107) 



This fixes the constant k = — 4/(3pp). If the pre- 
radiation era is made of BECs, this relation determines 
the ratio between the coupling constant A < of the 
bosons and their mass m (see Appendix IB") . 

When a > ai, we can make the approximation pj pp ~ 
(ai/a) 4 . This correspond to the radiation era described 
by an equation of state p = pc 2 /3. The conservation of 
Prada 4 implies that p P a\ = p r ad,oal. Writing p rad ,o = 
Qrad, oPo, where po is the present density of the universe, 
and using Eq. (47), we obtain 



1/4 



^H t P = 1.9710 



-32 



Using Eq. (52), we get 

/on \ V 4 



1 



1.6110 



oil 



V 8tt y x/Wp 
The value of the characteristic length a\ is 
oi = 2.61 10~ 6 m. 



(108) 



(109) 



(110) 



It provides a lengthscale that is intermediate between 
the Planck length lp = 1.62 10 -35 m and the present size 
of the universe ao = 1.32 10 26 m. It corresponds to the 
typical size of the universe at the end of the inflationary 
phase (or at the beginning of the radiation era). On the 
other hand, in the radiation era, Eq. t|105[) reduces to 
pc 2 ~ ppc 2 (T /T*) 4 , so we can identify ppc 2 /T 4 with the 
Stefan-Boltzmann constant a. Using Eq. (|4U) . this yields 

1/4 



■3 



T P . 



(Ill) 



Numerically, T* = l.llTp = 1.5710 32 K. 

Regrouping these results, we obtain the basic equations 
of the non-singular model 



1 



V = o/?(l - 4p/pp)c 



(112) 
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W = -(l-ip/p P ), q = l-2p/p P , (113) 



pp 



(114) 



the value a± = 2.6110~ 6 m. Using the approximate ex- 
pression (|119p of the scale factor during the inflationary 
phase, we obtain 



8tt 



1/2 



In 



t P . 



(120) 



Pp J \Pp 



1^ 



1/4 



(a/«i) 7 



[(a/a^ + ir 



(115) 



(116) 



The entropy (|68|) is given by 



4 / 30^^- N 

3 V 8tt 



7T 

15 

87 



1/4 



(^otp) 3 / 2 



(117) 



Numerically, S/ks = 5.04 10 87 . Finally, we note that 
K= (8 7 r/3) 1/2 ip 1 . 

As the universe expands from a — to a = +oo, the 
density decreases from pp to 0, the equation of state 
parameter if increases from —1 to 1/3, the deceleration 
parameter q increases from —1 to 1, and the ratio (c s /c) 2 
goes from —7/3 to 1/3. The temperature starts at T = 0, 
reaches a maximum and decreases to zero. 



B. The inflationary phase 

During the inflationary phase (a <C a±), the density 
has an approximately constant value p ~ pp. The corre- 
sponding Hubble parameter is 



H 



'8ttG 



V 3 



-pp 



1/2 



1/2 



1 

tp" 



(118) 



Numerically, t x = 23.2t P = 1.2510~ 42 s. This time cor- 
responds typically to the end of the inflation. Therefore, 
during the inflation, the radius increases exponentially 
rapidly from a, = lp = 1.62 10~ m at t = U = to 
a = ai = 2.61 10~ 6 m at t = h = 1.25 10~ 42 s. 

On the other hand, in the inflationary phase, the tem- 
perature is related to the density and to the scale factor 
by 



7>~G 2 ) i 1 Pp) 



7/4 



\ 1/4 / \ 7 

15 \ I a 



a i 



(121) 



Using the approximate expression (|119[) , we find that the 
temperature increases exponentially rapidly with time as 

T(t)~T P (§) V4 f^eW) 1 "'/". (122) 

When t < U = 0, the temperature is less than Tj = 
3.9110" 205 Tp = 5.5410 _173 K, and it tends to zero ex- 
ponentially rapidly when t — > — oo. Therefore, the pre- 
radiation era is extremely cold (this is consistent with 
the BEC model which assumes T = [13]). When 
t > U = 0, the temperature increases exponentially 
rapidly. It passes from T = T, = 5.5410- 173 K at 
t = U = to T = T x = 1.11 T P = 1.57 1 32 K at 
t = ti = 1.2510- 42 s. 

These values are based on the asymptotic results (|119[) 
and (|122p . They will be slightly revised in Sec. IVDI using 
the exact solution (|128[> of the Friedmann equations. 



Numerically, p = 5.16 10" g/m 3 and H = 5.37 lO^s -1 . 
Integrating Eq. (jTTSj) , we find that the scale factor in- 
creases exponentially rapidly with time as 



a(t) ~ Ipe 



(87r/3) 1/2 t/*i= 



(119) 



We have determined the constant of integration such that 
a(t = 0) = lp. This universe exists at any time in the 
past, so there is no singularity. 

When t < U = 0, the radius of the universe is smaller 
than the Planck length lp = 1.6210 _35 m and it tends 
to zero exponentially rapidly when t — > — oo. When 
t > U = 0, the radius of the universe increases exponen- 
tially rapidly on a timescale of the order of the Planck 
time tp = 5.39 10 _44 s. This phase of inflation connects 
the pre-radiation era to the radiation era on a very short 
interval of time (of the order of a few Planck times) . We 
call ti the time at which the scale factor a(t) reaches 



C. The radiation era 

After the inflation (a 3> ai), the universe enters in 
the radiation era. In that case, the quadratic pressure 
(oc — p 2 ) is negligible as compared to the pressure of ra- 
diation (p oc p) and we recover the standard radiation 
model of Sec. Ill E 21 The density is related to the scale 
factor by p ~ ppaf/a 4 . The Friedmann equation CO 
becomes 



H 



yielding 



lirG ppa\ ^ 1 



o \ i/ 2 1 - ,9 
S7T\ ' 1 /Ol\ 2 



tp v a 



32tt\ 1/4 /^ 1/2 



(123) 
(124) 
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Using the approximate expression (|124[) of the scale fac- 
tor during the radiation era, we find that the time at 
which a = a\ is 



32-k 



1/2 



t P . 



(125) 



Numerically, t x = 0.173t P = 9.3110" 45 s. This marks 
the beginning of the radiation era. Again, this numerical 
value will be revised in Sec. IV D[ We also have 



Pp 



3pp ( tp 



(a/ai) 4 



32tt V t 



and 



1 r \ 1/4 

15 \ ai 



T~T P — — ~ T, 



45 
32tt= 



1/4 



(126) 



During the radiation era, the density and the tempera- 
ture decrease algebraically as the universe expands. 



D. The general solution 



The equation of state (|112[) interpolates smoothly be- 
tween the pre-radiation era described by a constant den- 
sity p = pp (vacuum energy) and the radiation era de- 
scribed by a density p cx a~ 4 . It provides therefore a 
unified description of the early universe. Using the re- 
sults of Sec. IIV Dl the general solution of the Fricdmann 
equation (Tl4|) is f-47f : 

1 + v/(a/ai) 4 + A 



(a/ai) 2 
8^\ 1/2 i 



V(a/ai) 4 + 1 - In 



= 2 {t) r P +c > ^ 

where C is a constant of integration. It is determined 
such that a = lp at t = 0. Setting 



lp/ai, 



we get 



C(e) = y/^ + l - In 



(129) 



(130) 
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FIG. 9. Evolution of the scale factor a with the time t in 
logarithmic scales. This figure clearly shows the phase of in- 
flation connecting the pre-radiation era to the radiation era. 
During the inflation, the scale factor increases by 29 orders of 
magnitude in less than 10 -42 s. 
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FIG. 10. Evolution of the scale factor a with the time t in 
linear scales. The dashed line corresponds to the standard 
radiation model of Sec. Ill E 21 exhibiting a finite time singu- 
larity at t — (Big Bang). When quantum mechanics is taken 
into account (as in our simple model), the initial singularity 
is smoothed-out and the scale factor at t — is equal to the 
Planck length lp = 1.62 10 -35 m. This is similar to a second 
order phase transition where the Planck constant plays the 
role of finite size effects (see Sec. IVIII for a development of 
this analogy). 



Numerically, e = 6.20 1CT 30 and C = -134. For t 
— oo, we have the exact asymptotic result 



a(t) ~ ai e^^ 1/2t ^+ D , 



(131) 



with D = (C + In 2 — l)/2. Due to the smallncss of e, 
a very good approximation of C and D is given by C ~ 
1 — In 2 + 21ne and D ~ hie. With this approximation, 
Eq. HBP) returns Eq. (TiT9"j) , 
At t = U = 0, 



Pp 



(132) 



T t =T F 



1/4 



(e 4 + l)2- 



(133) 



Pi — Pp = 
5.5410" 173 K. 



Numerically, a* = l P = 1.62 1 - 35 1 
5.1610"g/m 3 , and T t = 3.91 10- 205 T P 

Using the results of Sec. IIV AI we find that the universe 
is accelerating for a < a c (i.e. p > p c ) and decelerating 
for a > a c (i.e. p < p c ) where 



1 



ai, 



Pc 



;Pp- 



(134) 
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FIG. 11. Evolution of the density p with the time t in loga- 
rithmic scales. During the inflation, the density remains ap- 
proximately constant with the Planck value p max = pp which 
represents an upper bound. 



FIG. 13. Evolution of the temperature T with the time t in 
logarithmic scales. For t < 0, the universe is extremely cold 
(T < 1CP 173 K). During the inflation, the temperature in- 
creases by 204 orders of magnitude in less than 1CP 42 s. Dur- 
ing the radiation era, the temperature decreases algebraically. 




FIG. 12. Evolution of the density p with the time t in linear 
scales. The dashed line corresponds to the standard radiation 
model of Sec. Ill E 21 exhibiting a finite time singularity at 
t = 0. Quantum mechanics limits the rise of the density to 
the Planck value p P = 5.16 10 99 g/m 3 . 



The corresponding temperature is 

1/4 



(135) 



The time t c at which the universe starts decelerating is 
given by 



1/2 



l r 



\/2 - ln(l + V2) - C] t P . (136) 



This is the time at which the curve a(t) presents an in- 
flexion point. It turns out that this inflexion point coin- 
cides with a\ (this property is true only for a = 1/3). 
Therefore, it coincides with our definition of the end 
of the inflation: t r = t\. Numerically, p c = p\ 



0.5 p P = 2.58 10" g/m 3 , a c = a x = 
Tj = 0.278 T P = 3.93 10 31 K, and t, 
1.2510" 42 s. 
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FIG. 14. Evolution of the temperature T with the time t in 
linear scales. The dashed line corresponds to the standard 
radiation model of Sec. Ill E 21 exhibiting a finite time singu- 
larity at t = 0. Quantum mechanics limits the rise of the 
temperature to about half the Planck value Tp = 1.42 10 32 K. 
In the pre-radiation era, the temperature drops to zero expo- 
nentially rapidly. 



The temperature starts from T = at t — > -co, in- 
creases exponentially rapidly during the inflation, reaches 
a maximum value, and decreases algebraically during the 
radiation era. Using the results of Sec. IIIIB1 we find that 
the point corresponding to the maximum temperature is 



7 1 /V, 



T. = I 



7/4 



15 

8^ 



1/4 



T, 



p- 



= It is achieved at a time 



U = 



8tt 



v/8-In 



l + vr 

V7 , 



c 



tp- 



(137) 



(138) 



(139) 
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Numerically, p e = 0.125 p P = 6.44 10 98 g/m 3 , a e = 
1.63ai = 4.2410~ 6 m, T e = 0.523 T P = 7.4010 31 K, and 
t e = 23.6 t P = 1.2710~ 42 s. 

The evolution of the scale factor, density and temper- 
ature as a function of time are represented in Figs. [911141 
in logarithmic and linear scales. 



E. Summary 

In this model, the universe starts at t — > — oo with 
a vanishing radius a — 0, a finite density p = pp = 
5.16 10" g/m 3 , and a vanishing temperature T = 0. 
The universe exists at any time in the past and does 
not present any singularity. It has been called "aion- 



iotic" universe in [47 1 since it has no origin nor end. 
To make contact with the Big Bang theory, we define 
the "original" time t — U = as the time at which 
the radius of the universe is equal to the Planck length. 
Thus di — lp — 1.6210 _35 m. The corresponding den- 
sity and temperature are pi ~ pp = 5.16 10" g/m 3 and 
% = 3.91 10~ 205 T P = 5.54 10~ 173 K. We note that quan- 
tum mechanics regularizes the finite time singularity that 
arises in the standard Big Bang theory. This is simi- 
lar to finite size effects in second order phase transitions 
(see Sec. IVII|) . We also note that the universe is very 
cold at t = ti = 0, unlike what is predicted by the Big 
Bang theory (a naive extrapolation of the law T oc i" 1 / 2 
leads to T(0) = +00). The universe first undergoes a 
phase of inflation during which its radius and tempera- 
ture increase exponentially rapidly [see Eqs. (|119[) and 
(|122[) ] while its density remains approximately constant 
[see Eq. (fTU7|) ]. The inflation "starts" at U = and 
ends at t\ = 23.3ip = 1.2510~ 42 s. During this very 
short lapse of time, the radius of the universe grows from 
a t = l P = 1.6210~ 35 mto 01 = 2.61 10~ 6 m, and the tem- 
perature grows from T t = 3.91 10~ 205 T P = 5.54 10" 173 K 
to Ti = 0.278 T P = 3.9310 31 K. By contrast, the den- 
sity does not change significatively: It goes from pi ~ 
p P = 5.16 10" g/m 3 to p x = 0.5 p P = 2.58 10" g/m 3 . 
After the inflation, the universe enters in the radiation 
era and, from that point, we recover the standard model. 
The radius increases algebraically [see Eq. (|124p ] while 
the density and the temperature decrease algebraically 
[see Eqs. (|126[1 and (|127[) ]. The temperature achieves 
its maximum value T e = 0.523 T P = 7.4010 31 K at 
t = t e = 23.6 t P = 1.2710~ 42 s. At that moment, the 
density of the universe is p e = 0.125 pp — 6.44 10 98 g/m 3 
and its radius a e = 1.63 ai = 4.2410~ 6 m. During the 
inflation, the universe is accelerating and during the ra- 
diation era it is decelerating. The transition takes place 
at a time t c = ti = 23.3 tp = 1.25 10 -42 s coinciding with 
the end of the inflation (i.e. a c = ax). 

We note that the inflationary process described above 
is relatively different from the usual inflationary scenario 
[37l [38| . In standard inflation, the universe is radiation 
dominated up to ti = 10~ 35 s, but expands exponen- 
tially by a factor 10 30 in the interval ti < t < tf with 



tf = 10 s. For t > tf the evolution is again radiation 
dominated. At t = U the temperature is about 10 27 K 
(this corresponds to the epoch at which most "grand 
unified theories" have a significant influence on the evo- 
lution of the universe) . During the exponential inflation, 
the temperature drops drastically; however, the matter 
is expected to be reheated to the initial temperature of 
10 27 K by various high energy processes Jfjij]. In the in- 
flationary process described previously, the evolution of 
the temperature is different. 

Remark: Of course, our definition of the "original" 
time t = 0, corresponding to a(0) = lp, is relatively 
arbitrary. What essentially matters is the evolution of 
the different variables on a tiny interval of time of the 
order of 20tp. 



VI. A MODEL OF NON-INFLATIONARY 
UNIVERSE WITH A NEW TYPE OF INITIAL 
SINGULARITY 



The equation of state (|55j) with n > and k > leads 
to a model of non-inflationary universe with a new type 
of initial singularity. We take the same values of a, n, 
|0*, ai and T* as in the previous section. We also take 
k = 4/(3pp). The basic equations of the singular model 
are 



P= ^(1 + 4p/pp)c 2 , 



w = -(l + 4p/p P ), 



1 + 2p/p P , 



Pp 



(a/ai) 



T = T f 
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1/4 



1' 



PP) \PP 



1/4 



T = T f 



1/4 



(a/ai) 7 



[(a/ ai ) 4 -ir 



(140) 



(141) 



(142) 



(143) 



(144) 



The entropy is still given by Eq. (|117j) . As the universe 
expands from a = ai to +00, the density and the temper- 
ature decreases from +00 to 0, the parameter w decreases 
from +00 to 1/3, the deceleration parameter q decreases 
from +00 to 1, and the ratio (c s /c) 2 decreases from +00 
to 1/3. 

For the equation of state (|140[) . using the results of Sec. 
IIV CI the general solution of the Fricdmann equation (TI"4"|) 
is EE SI: 



\/ (a/ai) 4 — 1 — arctan sj (a/a\) i — 1 



8^\ 1/2 t 



t P ' 

(145) 
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FIG. 15. Evolution of the scale factor a with the time t in 
logarithmic scales. There is no inflation since the universe 
already starts with a "large" radius ai = 2.61 1CF 6 m. How- 
ever, the initial condition is singular since the density and the 
temperature are infinite. 



FIG. 17. Evolution of the density p with the time t in loga- 
rithmic scales. The slope goes from —2/3 in the pre-radiation 
era to —2 in the radiation era. 




scales. The dashed line corresponds to the standard radiation 
FIG. 16. Evolution of the scale factor a with the time t in model of Sec. HIE 21 
linear scales. The dashed line corresponds to the standard 
radiation model of Sec. Ill E 21 (Big Bang). 



The universe starts from a singularity at t = 0. The 
initial radius of the universe is finite (with the value 
a(0) = a\ = 2.61 10~ 6 m) but the initial density and the 
initial temperature are infinite. The universe is always 
decelerating. The density and the temperature decrease 
indefinitely as the universe expands. There is no inflation 
since the universe starts with a "large" radius a\ ^> Ip. 
When t -> 0, 
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FIG. 19. Evolution of the temperature T with the time t 
in logarithmic scales. The slope goes from —4/3 in the pre- 
radiation era to —1/2 in the radiation era. 
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When t ^> tp, the universe enters in the radiation era. 
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the exponent 1/2 is the same as in mean field theories 
of second order phase transitions (i.e. M ~ (T — X^) 1 / 2 ) 
but this is essentially a coincidence. 

When quantum mechanics effects are taken into ac- 
count, as in the simple model of Sec. El the singularity 
at t = disappears and the curves a(t), p(t) and T(t) 
arc regularized. In particular, we find that a = l P > at 
t = 0, instead of a = 0, due to the finite value of Ti. This 
is similar to the regularization due to finite size effects 
(e.g. the system size L or the number of particles N) 
in ordinary phase transitions. In this sense, the classi- 
cal regime h — > is similar to the thermodynamic limit 
(L — > +oo or N — > +oo) in ordinary phase transitions. 



FIG. 20. Evolution of the temperature T with the time t in 
linear scales. The dashed line corresponds to the standard 
radiation model of Sec. Ill E 21 



In that case, we recover the results of Sec. IV CI 

The condition that the velocity of sound must be less 
than the speed of light imposes a maximum bound on 
the density. Using the results of Sees. IIII Dl and llV CI we 
find that the physical universe starts at a time 

/ 3 \ 1/2 1 

* i= V8W 2 [ 2 - arctan (2)]£p- ( 149 ) 
At this time, its radius, density and temperature are 

a,i = 5 1/4 ai, p l = \pp, (150) 

For t < ti, the solution (|145[) may be unphysical be- 
cause the velocity of sound is larger than the speed 
of light. Numerically, pi = 0.25 p P = 1.29 10" g/m 3 , 
c H = 1.50 oi = 3.9010~ 6 m, T 4 = 1.16 T P = 1.64 10 32 K 
and ti = 0.154 t P = 8.32 10" 45 s. 

The evolution of the scale factor, density and temper- 
ature as a function of time are represented in Figs. [T51l2"01 
in logarithmic and linear scales. 

VII. ANALOGY WITH PHASE TRANSITIONS 

The standard Big Bang theory is a classical theory in 
which quantum effects are neglected. In that case, it ex- 
hibits a finite time singularity: The radius is equal to 
zero at t = while the density and the temperature are 
infinite. For t < 0, the solution is not defined and we may 
take o = 0. For t > the radius increases as a a t 1 ! 2 . 
This is similar to a second order phase transition if we 
view the time t as the control parameter (e.g. the tem- 
perature T) and the scale factor a as the order parameter 
(e.g. the magnetization M). It is amusing to note that 
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FIG. 21. Effect of quantum mechanics (finite value of the 
Planck constant) on the regularization of the singular Big 
Bang solution (Ti — 0, dashed line) in the model of Sec. [Vj 
The singularity at t — is replaced by an inflationary expan- 
sion from the pre-radiation era to the radiation era. We can 
draw an analogy with second order phase transitions where 
the Planck constant plays the role of finite size effects (see the 
text for details). 




FIG. 22. Effect of quantum mechanics (finite value of the 
Planck constant) on the evolution of the scale factor a(i) in 
the model of Sec. I VII The standard Big Bang solution (h — 0) 
is represented by a dashed line. For Ti > 0, the initial radius 
is strictly positive and it tends to zero as h — > 0. 
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To study the convergence toward the classical Big Bang 
solution when Ti — > 0, it is convenient to use a proper nor- 
malization of the parameters. We call ho = 1.05 10 -34 J s 
the true value of the Planck constant (in this section, the 
subscript refers to true values of the parameters) . Then, 
we express time in terms of t P and lengths in terms of 
a?. We introduce a = a/a? and t = t/t P . Using lp = ctp 
and Eq. (|108[) , we obtain 

l P /l P =tp/t P = {a 1 /a° 1 ) 2 = h 1/2 , (152) 
where h — Ti/Tiq. Therefore 

e / £ o = 7?TT-o = a i/ a i = h ' ( 153 ) 
tp/a 1 

where eo = Ip/af = 6.20 1CP 30 . Using these relations 
and Eq. (|128p . we obtain 

/ ; / 1 + Jayh + A 

= 2 (f) 1/2 £A 1/2 + C(„» 1/4 ), (154) 

where C(e) is defined by Eq. ()130[) . This equation de- 
scribes a phase transition between the pre-radiation era 
and the radiation era. The normalized Planck constant 
h plays the role of finite size effects. Finite size scalings 
are explicit in Eq. (|154[) . For K = 1, we recover the 
non-singular model of Sec. [Vj For h = 0, we recover the 
singular radiation model of Sec. HIE 21 (Big Bang). The 
convergence towards this singular solution as h — > is 
shown in Fig. [2TJ 

We can make the same study for the model of Sec. IVII 
although the analogy with phase transitions with finite 
size effects is less clear since this model remains singular 
for finite values of the Planck constant. Nevertheless, we 
have represented the effect of quantum mechanics on the 
form of the curve a(t) in Fig. For finite values of h, 
using Eq. f|145[) , the evolution of the scale factor is given 
by 

yja 4 /h- 1 - arctan ^Ja 4 /h-l = 2 ^ i/h 1 ' 2 . 

(155) 

-1/4 

At t = 0, the scale factor is d(0) = Ti 

VIII. CONCLUSION 

In this paper, we have carried out an exhaustive study 
of the generalized equation of state (|55p for positive in- 
dices n > 0. For a = 1/3, n = 1 and k = —4/(3p P ), 
this equation of state describes in a unified manner the 
transition between the pre-radiation era (p = p P ) and 



the radiation era (p oc a -4 ) in the early universe. It pro- 
vides a model of early inflation. The case of negative 
indices n < is treated in Paper II. For a = 0, n = — 1 
and k = —p\, where p\ is the cosmological density, this 
equation of state describes in a unified manner the tran- 
sition between the matter era (p oc a -3 ) and the dark 
energy era (p = p\) in the late universe. It provides a 
model of late inflation. Combining these two approaches, 
we obtain a simple non-singular model for the whole evo- 
lution of the universe based on two symmetric polytropic 
equations of state (see Paper II). 

When we started our study, we were mainly interested 
in the equation of state p = (ap + kp 2 )c? introduced 
in our previous paper [471 ]. This equation of state was 
obtained in the context of Bose-Einstein condensates, 
assuming that they are the constituents of dark mat- 
ter. Although BECs should form after the radiation era 
when the universe has cooled sufficiently, we wondered 
whether this equation of state could have some sense be- 
fore the radiation era. In that picture, the pre-radiation 
era would correspond to primordial BECs with attractive 
self-interaction A ~ — (m/Mp) (see Appendix |B|) . Since 
the pre-radiation era is extremely cold (see Sec. EJ), the 
BEC model (which assumes T = 0) may be an interesting 
suggestion to develop further. 

However, we progressively realized that the study of 
a generalized equation of state of the form p = (ap + 
kp 1+1 ' n )c 2 was of interest beyond the context of BECs. 
Therefore, our present point of view is more general. We 
believe that the equation of state (f5"5")) may have sev- 
eral origins that should be discussed in future works. 
By studying this equation of state in full generality, we 
realized that the positive indices n > describe the 
early universe while the negative indices n < describe 
the late universe. Furthermore, a positive polytropic 
pressure (k > 0) leads to past or future singularities 
(or peculiarities) while a negative polytropic pressure 
(k < 0) leads to non-singular models. They exhibit 
phases of early and late inflation associated with a maxi- 
mum density p ma x = Pp (Planck density) corresponding 
to the vacuum energy in the past and a minimum density 
Pmin = PA (cosmological density) corresponding to the 
dark energy in the future. Therefore, in the polytropic 
model, the description of the early and late universe ap- 
pears to be very "symmetric" . This result is obtained 
in a purely theoretical manner, without reference to ob- 
servations. Strikingly, this symmetry is consistent with 
what we know about the real universe. This point will 
be further developed in Paper II. 

Another motivation of our study was to discuss some 
aspects of the history of science and publicize (in the 
same vein as [3, Q) the important contributions of 
Lemaitre in cosmology, which may not be sufficiently 
recognized. They include: 1. The first correct under- 
standing of the mysterious de Sitter model [l4j|; 2. The 
re-discovery of the Fricdmann equations and their ther- 
modynamical interpretation [151 : 3. The instability of 
the Einstein static universe [15l |: 4. The interpretation 
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of the redshifts as a results of the expansion of the uni- 
verse in accord with the Einstein theory of relativity fl5| ; 
5. The discovery of the Hubble law and the calculation 
of the Hubble constant two years before Hubble 6. 
The basis of Newtonian cosmology [2!| [58| ; 7. The cor- 
rect place of the cosmological constant in the Poisson 
equation Q; 8. The importance of the cosmological con- 
stant and its interpretation in terms of a vacuum energy 
density [58[; 9. The notion of a primordial singularity 
(primeval atom) that finally led to the Big Bang theory 
[27l I29I [30j ; 10. The imp ortance of quantum mechanics 
in the early universe [361 ] . Some other important contri- 
butions of Lemaitre have been collected recently in [65| . 



Appendix A: Equation of state p = (— p + kp 1 )^ 2 with 

n > and k > 
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FIG. 23. Evolution of the density and temperature as a func- 
tion of the scale factor a. We have taken n = 1. 



In this Appendix, we specifically study the equation of 
state (|55p with a = — 1, n > 0, and k > 0, namely 



p = (—p + kp J )c 2 . 



(Al) 



It generalizes the equation of state p = — pc 2 of the vac- 
uum. For the equation of state (|A1[) , the continuity equa- 
tion (fT2]) can be integrated into 



P* 



ln(a/a*) r 



(A2) 



where p* = (n/3k) n and a* is a constant of integration. 
The density is defined for a > a*. When a — > a*, p — > 
+00 and p — > +00; when a — > +00, p — > and p — > 0. 

The thermodynamical equation (|16|) can be integrated 
into 



T = T* 



(n+l)/n 



.Stp./P) 1 ^ 



(A3) 



where T* is a constant of integration. Combined with 
Eq. (|X2|) . we obtain 



ln(a/a»)" +1 \a 



(A4) 



When a — > a* and a — > +00, T — > +00. The temperature 
has a minimum at 



Pe 
P* 



n + l 



I 



n + 1 



Oe 

a * 



n+l 



= («+l)/3 



n+l 



(A5) 



(A6) 



The evolution of the density and temperature as a 
function of the scale factor is plotted in Fig. [531 

The equation of state can be written as p = wpc 2 with 



n o 
w = -! + - — 



3 \p* 



l/n 



(A7) 
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FIG. 24. Evolution of w, q, and (c s /c) 2 as a function of the 
scale factor a. We have taken n = 1. 



The pressure vanishes (u> = 0) at 

P ™/3 



P* 



On 

a* 



— = e" ( - 



n+l 



(A8) 

When a — > a*, w — > +00; when a — > +00, w — > — 1. The 
pressure is positive for a* < a < a w and negative for 
a> a w . 

The deceleration parameter is given by Eqs. ([77]) and 
([78]) . Together with Eq. (fAT)) . we obtain 



9= -1 



n /_p 
2 Vp* 



l/n 



(A9) 



The curve a(t) presents an inflexion point (a = q = 0) at 

n+l 



/9* 



a* 



Tc = (2 



= n/2 ±c _ f_ \ 3n/2 



(A10) 

When a — >• a*, 5 — > +00; when a — > +00, g — > — 1. The 
universe is decelerating when a* < a < a c and accelerat- 
ing when a > a c . 
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Finally, the velocity of sound is given by 

1/n 



-§ = -! 



1 



(All) 



The velocity of sound vanishes at the point defined by 
Eqs. (|A5I) - (|A6|) at which the temperature reaches its 
minimum. At that point, the pressure is minimum with 
value 

At = "7 — ^— rr • (A12) 

When a—t a*, (c s /c) 2 —> +oo; when a — » +oo, (c s /c) 2 —> 
— 1. The velocity of sound is real when a* < a < a e and 
imaginary when a > a e . On the other hand, the velocity 
of sound is equal to the speed of light at 
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Zi = P ("+l)/6 



(A13) 
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FIG. 25. Evolution of the scale factor with time. We have 
taken n = 1. For this index, KU = 0.128, Kt c = 0.236, and 
Kt c = 0.363. 
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For n = 1, this corresponds to the point at which the 
pressure vanishes (w = 0). The velocity of sound is larger 
than the speed of light when a* < a < a s and smaller 
when a > a s . The evolution of w, g, and (c s /c) 2 as a 
function of the scale factor a is represented in Fig. |2"41 

Setting i? = a/a*, the Friedmann equation (|14[) can be 
written 

R = KR ,„ , (A15) 



where K 



(lni?)«/2' 
•SttG/0,,,/3) 1 / 2 . Its solution is 

^+2 ^^ 2 /(«+2) 



(A16) 



where we have determined the constant of integration 
such that R = 1 at t — 0. The density and the tempera- 
ture evolve as 

-2n/(n+2) 



-Kt 



-Kt 



-2(n+l)/(n+2) 



(A17) 



.(A18) 



The universe starts at t = with a radius i? = 1 and an 
infinite density and infinite pressure. Then, the universe 
expands indefinitely while the density decreases to zero. 
The physical solution, for which the velocity of sound is 
smaller than the speed of light, starts at t = t{ where 
Kti = [2/(n + 2)][(n + l)/&\ {n+2 ^ 2 . For n = 1, this 
corresponds to the time at which the pressure becomes 
negative. The expansion of the universe is decelerating 
when t < t c and accelerating when t > t c , where Kt c = 
[2/(n + 2)](n/2)("+ 2 )/ 2 . The velocity of sound becomes 
imaginary when t > t e , where Kt e = [2/(n + 2)][(n + 
l)/3](™+ 2 )/ 2 . This corresponds to the time at which the 
temperature reaches its minimum. The evolution of the 
scale factor with time is plotted in Fig. [231 



Appendix B: Connection to the BEC model 

At T = 0, and in the strong coupling limit, the equa- 
tion of state of Bose-Einstein condensates (BECs) with 
short-range interactions is 



2ira*h 



P 



-P 



(Bl) 



where m is the mass of the bosons and a s is the s- 
scattering length (43|. Therefore, a self- interacting BEC 
has a polytropic equation of state p = Kp 1 with in- 
dex 7 = 2 (i.e. n = l) and polytropic constant 
K = 2Tra s h 2 /to 3 . Furthermore, the scattering length 
a s can be positive or negative (positive scattering cor- 
responds to repulsive self-interaction, and negative scat- 
tering to attractive self-interaction). In terrestrial BEC 
experiments, some atoms like 7 Li have a negative scat- 
tering length [6o| . As a result, the pressure of a BEC 
can be negative. Therefore, there exist physical systems 
in which it is possible to obtain polytropic equations of 
state with negative pressure. Furthermore, the origin of 
these equations of state is due to quantum mechanics 
since the Planck constant explicitly appears in Eq. (|B1[) . 
This may be a justification (but not the only one) for the 
generalized equation of state ([55)) . 

It is convenient to define a coupling constant A by [39| : 



A a s mc 
8tt ~ h 



(B2) 



We can now use (to, A) as independent variables instead 
of (to, cl s ) . The equation of state (|Bip can be rewritten 

as 



V 



j 2 2 

hp c , 



with k 



Am 4 c 3 ' 



(B3) 



In the inflationary model of Sec. |V1 in which the pre- 
radiation era is described by a "fluid" with a polytropic 
equation of state of the form (|B3I a) with k < 0, the 
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maximum density is p* = 4/(3 \k\). If we assume that 
this "fluid" is made of BECs, then, using Eq. (|B3I b). we 
obtain 

16 / m \ 4 1 

p* = Kpp, With K = — I — — j—r. (B4) 



Assuming that = p max is of the order of the Planck 
density pp, i.e. k ~ 1, we find that the coupling constant 

of the BEC scales as A (16/3)(m/M P ) 4 . Of course, 

the nature of these "primordial" BECs remains highly 
speculative. It is possible that this BEC description is 
just an effective one. 
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